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Abstract 

We show that supersymmetric "new physics" beyond the MSSM can naturally accom- 
modate a Higgs mass near 126 GeV and enhance the signal rate in the h — > 77 channel, 
while the signal rates in all the other Higgs decay channels coincide with Standard Model 
expectations, except possibly the h — > Z7 channel. The "new physics" that corrects the 
relevant Higgs couplings can be captured by two supersymmetric effective operators. We 
provide a simple example of an underlying model in which these operators are simultane- 
ously generated. The scale of "new physics" that generates these operators can be around 
5 TeV or larger, and outside the reach of the LHC. 
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1 Introduction 



The ATLAS and CMS collaborations at CERN recently presented strong experimental ev- 
idence for a Higgs-like resonance around 126 GeV [1], marking a historic achievement in 
particle physics. The signal rates in the ZZ* and WW* channels are in good agreement with 
the Standard Model (SM) predictions. The bb and r + r~ signal rates are also compatible 
with the SM, although with substantial error bars. At the time of writing, the signal rate in 
the h — )■ 77 channel is about 1.5 — 2 times larger than the SM prediction. This discrepancy 
is only at the level of two standard deviations and there are theoretical uncertainties [2]. 
Nevertheless, it is interesting to contemplate whether physics beyond the SM (BSM) can be 
responsible for this discrepancy (for recent work in this direction, see [3]). 

In this work we shall assume that the excess in the diphoton channel is due to BSM 
physics that has negligible effect on channels other than h — > 77. Indeed, this channel is 
sensitive to new physics, since it is a loop-level process in the SM. With this in mind, we 
shall focus on the minimal supersymmetric (SUSY) extension of the SM (MSSM) close to 
the "decoupling limit" [4], in which the lightest neutral CP-even Higgs boson h is SM-like. 

It is surprising that, despite its large number of parameters, the MSSM has difficulties in 
accommodating an enhancement of the h — > 77 partial decay width I\ 77 without affecting 
the other partial decay widths. In fact this requirement seems to single out loop- induced 
contributions from very light color singlet superpartners with a significant coupling to the 
Higgs, meaning strongly mixed light stau sleptons, at around 100 GeV [S]. However, this 
introduces issues with vacuum stability and may even be possible to rule out at the LHC. 
In addition, large radiative corrections are needed to obtain a mass of ~ 126 GeV for 
the lightest Higgs of the MSSM. This requires large supersymmetry breaking terms, such as 
TeV stop masses and/or a large top A-term. The lack of evidence for superpartners in the 
direct SUSY searchers at the LHC also indicates that soft terms should be large. However, 
large supersymmetry breaking terms lead to severe fine-tuning [61 [7] in most versions of the 
MSSMj^] (with or without universal gaugino masses, or Higgs soft masses different from mo 
and also from each other). This situation suggests that a solution to this problem is not in 
the SUSY breaking sector but rather in the one that preserves supersymmetry. Following this 
idea, the purpose of this work is to answer whether one can have a minimal supersymmetric 
extension of the Higgs sector that allows ~ 126 GeV without undue fine tuning and a 
simultaneous h — )■ 77 enhancement, while complying with the negative SUSY searches so far. 

To this end, one way to proceed is suggested by minimal extensions of the MSSM Higgs 
sector, like the Next-to-MSSM (NMSSM) model which contains an additional singlet chiral 
superfield (see [9] for a review), and by models where the soft B[i term is promoted to a 
SUSY operator [10]. In the NMSSM an enhancement of the h — > 77 branching ratio is 
possible (although this also alters other couplings) beyond the SM level [11]. However, the 
NMSSM remains badly fine-tuned (fine tuning A > 200) for m h « 126 GeV p2]. There are 

1 For a discussion of the negative impact of the EW fine tuning on the \ 2 fit of such models see also [8]. 
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known ways to bypass this problem such as in the so-called "generalized" version of NMSSM 
(GNMSSM) with a superpotential mass term for the singlet superfield, where the electroweak 
fine tuning is significantly reduced to more acceptable levels (A rj 30) [121 E31 El • These are 
examples of explicit supersymmetry-preserving modifications of the MSSM that can render 
more natural the interpretation of the resonance at 126 GeV as the lightest Higgs. 

In this work, instead of considering such specific models, we shall relax the rigid, minimal 
structure of the MSSM Higgs sector and perform an effective field theory analysis of the 
most relevant SUSY-preserving operators in this sector. This approach should recover, in a 
particular region of the parameter space, scenarios such as those presented above (GNMSSM, 
etc). We show that it is possible to naturally accommodate a Higgs mass of 126 GeV, an 
enhanced Higgs coupling to photons, and simultaneously SM-like Higgs couplings to the other 
particles, using only a few supersymmetry-preserving operators with small coefficients. 

In general there is a large set of operators that one could consider in the Higgs sector 
|15j . Regarding the Higgs mass it is known that the presence of the effective dimension 5 
superpotential operator [161 E! 

^{H u -H d f (1.1) 

can accommodate ~ 126 GeV without undue fine tuning [13] . The suppression scale M 
represents the mass scale of the SUSY degrees of freedom that have been integrated out to 
generate (1.1). Concerning T^-y-y, from the list of effective operators of dimensions d = 5 and 



d = 6 in the Higgs sector [131 02] one notices the presence of a SUSY effective operator 

^(H u .H d )Tr(W a W a ) (1.2) 

that can significantly modify the h — > 77 rate. Here W a is the electroweak gauge field 
strength superfield. Based on these observations we intend to investigate closer the phe- 



nomenological impact of (1.1) and (1.2). As a simple example will show, these operators can 
be generated simultaneously by an underlying microscopic model, making their combination 
rather natural. The effect of both operators is maximized for small tan /3, where the MSSM 
tree level contribution to the Higgs mass is minimized. This means that the impact of the 
d = 5 operator in accommodating an ~ 126 GeV is rather significant. To our knowledge, 



the particular combination of the effective SUSY operators in (1.1) and (1.2) has not been 
studied in the past for the problems we address 

The paper is organized as follows. In Section [2] we calculate the corrections from the 
effective operators to the Higgs mass and mixing angle. In Section [3] we discuss how these 
operators correct the Higgs couplings and signal rates, with focus on the decoupling limit. 
The results in terms of the Higgs mass and the partial widths for the h — > 77 and h — > Z7 
channels are discussed in Section [4} Section [5] provides an example of the origin of the effective 



The operator in (1.2 1 was separately discussed in [19] 
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operators and Section [6] contains our conclusions, while some details concerning the on-shell 
Lagrangian are given in the Appendix. 



2 Corrections from SUSY operators to the Higgs couplings 

The effective model we consider consists of the usual MSSM Higgs sector, extended by the 
operators discussed in the introduction. The relevant part of the Lagrangian is, in standard 
notation, 

C = Jd A 9 (1 - m 2 9 2 9 2 ) Hje v 'Hi + (Jd 2 9fi(l + B9 2 ) H d -H u + h.c.) + 5 + O e (2.1) 



i=u,d 



where the chiral superfields have components Hi = (hi,ipi, Fj), and mi and B are the soft 
terms. O5 is the only operator of dimension 5 that one can write in the Higgs sector, up to 
non-linear field redefinitions [20], and has the form: 



5 M 



f d 2 6 (H u ■ H d ) 2 + h.c. (2.2) 



For the component fields expression of O5 see eq. (A.3). 

There is a long list of operators in the Higgs sector of dimension d = 6 [1 3|. [T8l [20] . A 
careful analysis of these operators shows that of all these there is one of them that can couple, 
in a supersymmetric way, to two gauge bosons: 

s=l,2 "ssj 

Here g\ and 32 denote the U(l)y and SU(2)^ gauge couplings, respectively, and k s is a 
constant that cancels the trace factor. W a is the SUSY field strength of the U(l)y (SU(2)i) 
vector superfield V\ (V2) of components (A s , V S;/X , D s /2), s = 1,2. O5 and Oq provide a 
minimal set of operators that is enough for our purposes. One can also consider SUSY 



breaking effects associated to these operators (see the appendix, eq. (A.3)), but we only 
seek supersymmetric solutions to our problem. The effective expansion is reliable when 
c o,i,2 = 0(1) an d M is the largest scale in the theory. One can choose one of 00,1,25 for 
example Co, and set it to c$ = 1 by redefining M. But it is useful to keep Co to easily trace 
or turn off the effects of O5. Also, to modify the diphoton rate c\ or C2 (or a combination 
thereof) is enough, so together with the scale M we effectively have only two parameters. 
Keeping both c\^% generates an additional interesting coupling, see later. 

Additional operators of d = 6 can be present. Although they could have an impact on 
the Higgs mass |18j , they have an additional scale suppression relative tcP] O5 . There is an 



3 Strictly speaking, this is true for small tan j3 region, that will actually be the relevant region in our case. 
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operator similar to Oq but involving instead the SU(3)c gauge group, that we do not consider 
here; this would change dramatically the Higgs decay rate to gluons, away from the SM values. 
We take the good agreement with the SM in most channels as evidence that if present, the 
coefficient of this operator must be small. Finally, another reason to restrict our analysis to 
0§£ is that, as discussed later, they can be simultaneously generated by underlying physics. 



2.1 The on-shell Lagrangian 

The calculation of the on-shell Higgs Lagrangian extended by C5 and Oq is detailed in the 
appendix. The result is 



C 



D ' m (i + ^C^-^ + h-c.)) +(2-> I] 



CO 



H + 2 — h d ■ h u 



(N 2 + i^i 2 ) + [f {w x " 2 x > + w xl ) ^ hd? + |/l « |2 ) + kc -] 

+ { 7^2 • hd) [i {K^t^l - V^X a 2 )] + h.c. + (2^1)} 
+ ~M t 2 ( hu ' hd)($d ■ i>u) - (K -tpd + ipu- h d ) 2 ] +h.c. 



+ 



+ 



{ 



_C2_ 

AM 2 



\ (h u ■ h d ) {F^F^ V + l - e^F 2 a ^F 2 % a ) 



V2 (h u • + ^ • h d ) o» v \ a 2 F% av -^u-^d A2A2 ] + (2 -> 1) + h.c.} 



fiB (h d ■ h u ) + h.c. - fh 2 d \h d \ 2 - fh 2 u \h u \ 



(2.4) 



where fh 2 = m 2 + \fi\ 



u,d. For the explicit form of D 2 D 2 and Df, see eqs. (A. 7), (A. 8). 



Eq.(2.4) contains all the information one needs to extract the corrections to the Higgs 



masses and couplings. In particular, notice the presence of new, supersymmetric couplings: 



1 



(K • h d 



Tr Fi + 



ci 



TrFr 



A* + 2-^ h d ■ h u 



(\h d \ 2 + \h u \ 2 ) + h.c. (2.5) 



\M 2 2 M 2 

which are important below. There are also direct Higgs-higgsino and higgsino-gaugino cou- 



plings that can be relevant for dark matter models. From (2.4) we find the Higgs scalar 
potential Vh 

2 

V h = m 2 d \h d \ 2 + m 2 u \h u \ 2 - [fiBh d ■ h u + h.c] + \h\h u \ 



1 + 



1 



+ ^(IM 2 -I^l 2 ) 2 9 2 +[{h d -K) + ^ +h. c . 



+ 



M p 



gjCl g2 c 2 

M 2 M 2 
\h d \ 2 +\h u \ 2 )(h d -h u ) + h.c. 



+4 



Si 
M 

(9 2 



_C2_ 

2M 2 
2 



(h d ■ h u + h.c.) 
h d - h u \ 2 (\h d \ 2 + \h u \ 2 ) 



g 2 + g 2 2 ), (2.6) 



which depends on two parameters: cq from the effective dimension 5 operator and the com- 
bination (gfci + g 2 C2) from the effective dimension 6 operator. Note that last term in the 
first line above does not contribute to the neutral Higgs sector masses. 
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We also include dominant loop corrections, although they do not play the same crucial 
role they do in the MSSM. In the small tan/3 regime and for dominant top Yukawa coupling, 
the one- loop and leading two-loop correction to 14 is |21j . 

,2 



A 14 



a' 



5\h u \ 



(2.7) 



where 



3/4 
g 2 vr 2 



, Mr X t 1 

ln — + -r + 

7Tlt 4 



ox/ Mr\, Mr' 

v i o i n * ' In - 



16 5 3 X* + 21n 



2 (A t - //cot/3) 5 
M? 



(A f -/xcot/3) ; 
12 M| 



(2.f 



with M~ = m^m^, and 53 is the QCD coupling. 



2.2 The Higgs mass and mixing angle 



The scalars receive mass corrections from the usual one-loop radiative corrections but now 
also from the effective operators. Here we take the parameters co,ci,C2 to be real. We find 
the following result for the mass of the lightest Higgs scalar h: 



m 2 h 



- jm^ + m 2 z + 5 m\ sin 2 (3 — + Am| 



(2.9) 



where 



w 



( m A ~ m z) cos 2/3 + Sm z sin /3] + sin 2f3(m A + m z 



2 \2 



(2.10) 



and where Am^ is the contribution due to the higher-dimensional operators: 



Ami 



2/i 



£0 
M 



*\ + '■ 2 /' y ) S 2 + 



M 2 + M 2 



S3 + O 



\M 3 



(2.11) 



with 



si 



sin 



2/3 { 



1 + 



(m 2 ^ + m. 2 



(2.12) 



«2 



«3 



4// 2 



sin 2 2/3 + 



1 + 



2/i 5 



(ml + m|) sin 2 2/3) + — - (m\ + m|) 2 v 4 sin 2 2/3 



W 



7?^ sin 2/^ 

— sin 2/3 + [8m 2 , - (4 + 3<5)m| + 65m| cos 2/3 + 3(4m 2 4 - 5ml) cos 4/3] 

32 128 \/w 



where we kept (small) effects from the interplay between the effective operators and the 
one- loop correction to 14- The mass of the CP-odd Higgs boson is 



m\ 



2Bn 2v 2 /c 
sin 2/3 ~ sin 2/3 \M 



t> 4 cos 2 2/3 ( g{c\ g\c2 



32 sin 2/3 V M 2 M 2 



(2.13) 
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which does not receive one-loop corrections. The mixing angle a is given by 



with 



tan 2a 



V 

fz 
Sa 



1 

V 



2v' 



cq 



(mi+m|)tan2/3-^^ M 

/gig. g|c 2 \ t; 4 (m| /z -m\f A ) 
1 ' ' ' r ' " 1 32 P cos 2 2/3 



2^) t; 4 tan 2/3 
M/ 



M 2 



m\-m% + (sec 2/3 - l)<5m|/2 

4 cos 2(3 - (2 - 5 cos 2/3 + 6 cos 4/3-3 cos 6/3) 5/4 

cos 2/3 + 3 cos 6/3 



(2.14) 



(2.15) 



One can see the corrections to tan 2a due to the effective operators, that are used below. 

We also note that the new operators correct the gauge field kinetic terms when the Higgs 
fields receive vevs. The corrected gauge couplings are the ones that appear in the following. 

3 Corrections to the partial widths of h — > 77 and h — > Z7 

In this section we study how the new operators correct the Higgs couplings to the SM particles. 
To this end, we parametrize these corrections in terms of the usual MSSM Higgs couplings. 



3.1 Higgs couplings and signal rates 

The renormalizable part of the Lagrangian for the lightest neutral CP-even Higgs scalar h 
can be written 1221 as 



mt m c _ mb - 

£ren = ~ C t htt — C c hcC — Cb hbb 

V V V 

+c z m ^hZ^ Z„ + cw^h W + » W~ 
where the dimensionless coefficients are given by, 



c T — n T 
v 



cos a 



ct 



sin P 



sin a 



cos f3 



, cz = cw = sin(/3 — a) . 



(3.1) 



(3.2) 



where the mixing angle a is given in (2.14). In the scenario under consideration all loop 



corrections to the tree level coefficients in (3.2) are negligible. The usual SM values for 
the couplings in (3.1) and (3.2) are obtained in the decoupling limit, in which a — > /3 — 

1, where i = 



7r/2, implying that cos a — > sin/3, sin a 



cos /3 and hence, c. 



t,c,b,r, Z,W. 



We work in the limit when loop contributions from superpartners and other Higgs scalars 
are negligible. The dimension-five part of the Higgs Lagrangian, which takes into account 1- 
loop contributions from SM particles as well as the contributions from the effective operators 
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in (2.4) can be written as 

•£dim5 _ 



12-nv 

+ 1 < y 7 + < 



) 8irv 



flU 



JBMSSM^ "EM hz(1 „ Fi 



) 4-7T shift 



(3.3) 



The 1-loop contributions to these coefficients are given b}|^] 

>°p = ^Af + cfc^W » 1.03 q - (0.05 + 0.07i) c fe 



„loop 



-Z 7 



c^A^ + ct^W « -8.36 c w + 1-84 q 



5.80 cm/ - 0.31 ct 



(3.4) 



where, in the last steps, we have inserted 
explicit expressions are given in appendix 



A.2 



126 GeV in the 1-loop form factors A, whose 
In the decoupling limit, where Cj — > cf M = 1 
in (3.2), the c loop -coefficients in (3.4) approach the values they have in the SM, which, for 
rrih = 126 GeV, follow trivially from (3.4), 



„loop 

Joop 

"7 

- 7 Z 



SM 

SM 

SM 
c ~fZ 



0.98 + 0.07i 
-6.52 
5.49 . 



(3.5) 



In order to obtain the <J^ MSSM and c^ SSM coemcients in (2.2), we extract the following 



component interactions from the operators in (2.3) 



Ob D 



E ° s h .hAF a ^ v F a -\-- f ^P a F a F a l+hc 

s=l,2 V ~ / 



(3.6) 



-> V C °g^ 2 + a) ([ci cos 2 9 W + c 2 sin 2 9 w ]hF^F^ + 2(c 2 - ci) sin0 tt cos9 w hF^Z^ 



where we have used, 



0z,0 



h u -h d = hZh d - h u h 



K. 



o 



1 



>i + Re/i° + i!m/i°) , 



Re h d = — sin ah + cos a If , Re /t u = cos ah + sin a if 



Ai„ = cos 8 w An - sin w Zu 



A^l = sin W A u + cos6 w Z u 



(3.7) 



and t>d = wcos/3, v u = v sin/3, with u = 246 GeV. Moreover, the hypercharge gauge boson 

(3) 

Aifj, and the (third component of the) SU(2)^ gauge boson A 2 ^ have been rewritten in terms 
of the photon An and the Z boson Zn. Note that there is also a dimension 5 operator 



generated from (2.3) that involves the Higgs scalar h and two field strengths of the Z boson 



(as well as an analogous operator involving two field strengths of the W boson). However, 
since these operators will have couplings comparable to the 77 or Z7 couplings, but strongly 

4 See also [24] for additional studies of h — > Zj. 
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phase space suppressed, we expect them to be irrelevant with respect to the usual dimension 
3 Higgs coupling to the Z and W bosons in (3.1). Therefore, we do not consider them. 
The contributions to (3.3) from (3.6) are given by: 



^BMSSM 

"7 



BMSSM 
c -yZ 



7TV 2 COs(/3 + Q) 2 . 2 

(ci cos o w + C2 sm U w ) 



M 2 a E M 

7T v 2 cos(/3 + a) 



M 2 a E M 



(C2 - c\ ) sm 



, cos I 



(3.8) 



In the decoupling limit, where cos(/3 + a) 
maximized for small tan/?. 



sin 2/3, we see that the coefficients in (3.8) are 



We can now define the relevant Higgs partial decay widths, normalized to the correspond- 



ing SM value, in terms of the dimensionless c-coefficients in (3.2), (3.4), (3.5), (3.8) 

loop 



*■ hii 

pSM 

hii 



1 hgg 
' -pSM 
hgg 



„1°°P 



.SM 



1 /i77 

' pSM 

/i77 



+ c; 



BMSSM 

7 



.SM 



1 h-yZ 
L h-yZ 



loop , BMSSM 
c -yZ c -yZ 



.SM 

- 7 Z 



(3.9) 



as well as the corresponding branching ratios (BRs), 



BR 



/177 



nnSM 



BRfeii 
Biff 



„loop 



+ c 



Ctot 

BMSSM 

7 



C SM Ctot 



BR 



^99 



nnSM " 

hgg 

BRfe 7 z 

Tip SM 



.loop 



.SM 



Ctot 



loop , BMSSM 
C 7 Z C 7-^ 



C^Ctot 



(3.10) 



The coefficient c to t in (3.10) can be written as, 

ictoti 2 = Yl i 



BR^ + 



.loop 



.SM 



nnSM 
V&hgg 



(3.11) 



where we have neglected the contributions from, for example, h — )■ 77 and /i — > Z7, as well 
as possible invisible decays. Let us now define the inclusive, as well as the individual gluon- 
gluon fusion (ggF), vector boson fusion (VBF) and vector boson associated (VH) production 
cross sections, normalized with respect to the corresponding SM values, 



Q"incl 

°"ggF 
SM 
°ggF 



Joop 



,SM 
9 



T SM 
7 ggF 



+ \cy\ 2 (cr 



VBF 



+ rr SMN | 



a 



SM 



SM 
VBF 



+ 0" 



SM 
VH 



OVBF 

_SM^ 
VBF 



CVH 



SM 
VH 



(3.12) 



where we have denoted cy = cz = cw , since the Higgs couplings to Z and W bosons coincide 
in (3.2). We can now write, for example, the signal rates in the inclusive and dijet channels 
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of the h — )■ 77 decay mode, again normalized with respect to the SM, 



R 



incl 

'77 



ctot BRft. 



77 



(7 



SM 
tot 



BR 



SM 
A77 



odijct 
L 77 



J°°P /.SM 
'9 / c q 



' SM , 7 I |2 SM , 7 U^I^SM 

ggF ' e VBF 1^1 a VBF + e VH l c V| ^VH BR/ ( , 



77 



7 SM , 
e ggF°ggF + 6 



7 



SM 



VBF VBF 



SM 



VH^VH 



nnSM 



(3.13) 



where the e 7 -coefhcients are the selection efficiencies for the different production modes in 
the dijet-tag category of final states. 



3.2 The decoupling limit 

Let us now take the decoupling limit, in which, 



.SM 



„1°°P 



.SM 



„1°°P 



.SM 



^1 "7 
where i = t,c, 6, r, Z, W. This implies that |ctot| 



Joop 
J3M 



1 



1 in (3.11 ) and that, 



whereas, 



-L /177 

^SM~ 

/177 



77 



nnSM 



pSM 

hii 



1 + 



T3TJ SM 



.BMSSM 
"7, dec 



-pSM 
hgg 



R r>SM 



1 



„SM 



r SM 
1 h^Z 



BRfc 7 z 

nnSM 



1 + 



.BMSSM 

J 7Z,dcc 



.SM 

- 7 Z 



for which the coefficients in (3.8) are given by, in the decoupling limit, 

7r tP 1 sin 2/3 



BMSSM 
L 7,dec 



„BMSSM 
^7Z,dec 



M 2 a E M 

7tw 2 sin 2/3 

M 2 a E M 



(ci cos 2 ^ + C2 sin 2 
(02 — ci) sin 2 W cos ^ 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



In the decoupling limit, the production cross sections in (3.12) are all equal to their SM 
corresponding SM value, 

Omcl CTggF CVBF CVH 



a 



SM 



'incl ^ggF "VBF "VH 

Thus, all signal rates, for any production mode, associated with the channels h —> ii, for 
i = t,c, b, t, Z, W, as well as h — > gg, will be equal to their corresponding SM value. In the 



SM 



SM 



SM 



1 



(3.18) 



h — > 77 channel, we see that the signal rates in (3.13) (as well as any other signal rate in the 



h — > 77 channel) will be given by the corresponding normalized partial width, 



R 



77 



R 



incl 

77 



odijct 



1 /177 
T^SM" 

/177 



1 + 



^BMSSM 
"7, dec 



.SM 



(3.19) 
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and the same for the h — > Z7 channel, 



Tjincl 
K Z~f 



R 



.dijet 
Z 7 



r 



hZ-y 



1 hZy 



1 + 



„BMSSM 
^7, dec 



-Z7 



(3.20) 



In summary, in the decoupling limit, all the partial decay widths, except for r^ 77 and 
r^2 7 , are all equal to their corresponding SM value. This implies that all the production 
cross sections, as well as the signal rates in all other channels, are equal to their SM values. 
Moreover, as seen in (3.19) and (3.20), the partial decay widths for h — > 77 and h — > Zj, 



normalized with respect to the SM values, coincide with the corresponding signal rates. 
Hence, in this limit, Vf^/T^f and F 



hZ-y 



A^z 7 can be compared directly to the measured 



signal rates. From (3.17) notice that if c\=C2 one can change T] lll without affecting T^Zj- 



4 Results 



We can now evaluate the effect of the operators in (2.2) and (2.3) on the mass of the lightest 



neutral CP-even Higgs particle h and on the partial decay widths r^ 77 and T^z^, that directly 
correspond to the rates in the decoupling limit, as discussed in the previous section. 

The Higgs mass in (2.9) as a function of tan/3 is displayed in figure [T] It is well known 
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Figure 1: The mass of the lightest Higgs particle h as a function of tan/3, for M = 5 
TeV, ix = 300 GeV, m A = 1 TeV, no mixing (X t = 0), M~ t = 500 GeV (left panel) 
and Mj = 1 TeV (right panel), and for various values of the O5 coefficient cq. The 
solid curves include the MSSM loop corrections. 



that we can accommodate a Higgs mass at 126 GeV by tuning the soft parameters in the 



loop correction (2.7), but this usually demands a large tan/3, which we do not consider here 
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Figure 2: The h — )■ 77 signal rate i? 77 (solid black lines) in (3.19) and the h — > Z7 



rate Rz-y (dashed blue lines) in (3.20) are shown as functions of the coefficients c\ 
and C2 of the operators in ( |2.2[ ) and (2.3). In the plots we have set tan/3 = 3 (left 
panel) and tan/3 = 7 (right panel), M = 5 TeV and taken the decoupling limit. 



(since then additional Yukawa couplings that we do not include become important). With 
the d = 5 operator (2.2 ), one can easily obtain a value of mh ~ 126 GeV, see figure [T] this has 
an acceptable fine tuning A < 30 [13J , even for small tan /3 < 10, which is an otherwise very 
fine-tuned region of the MSSM. Therefore the effect of the d = 5 operator is more important 
than usually thought. 

The dimensionless parameter e = co/u/M measures the extent to which the contribution 
from 05 to the mass can be considered perturbative, and for the given numbers it is below 
e < 0.06. The mass contributions from Oq rarely affect noticeably the curves in figure [TJ but 
they are included for completeness (in the figure, c\ = C2 = — 1). We included only a subset 
of the loop corrections, that is relevant at low tan /3, so we expect the curves to differ from the 
complete result by few GeV only, which we confirmed in our examples using FeynHiggs |25j . 

For the value of the scale M found above in the region of 5 TeV, not within the LHC 
reach, one would like to examine the signal rates for h — > 77 and h — > Zj. These signal rates 
i? 77 and Rz-y in the decoupling limit, given in eqs. (3.19) and (3.20), are shown in figure [2] 



as functions of the coefficients c\ and C2 of the operators in (2.2) and (2.3). Concerning the 



h — > 77 channel, from the dependence on c\ and C2 in c^^. 81 ^ in (3.17), and from the fact that 



c 8M is negative in (|3.5|), we see that the maximal enhancement of i? 77 in (|3.19l) is obtained 



"7 



for negative values of both c\ and C2. In contrast, for positive (and not too large) values of 
the two coefficients, the h — > 77 signal is depleted with respect to the SM prediction, as can 
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Figure 3: The /i — )• 77 signal rate i? 77 from (3.19 ) as functions of the scale M, where 



we have set the coefficients of the two operators in (2.3) equal, c\ = C2. We have set 



Mi = 1 TeV, fi = 300 GeV, rriA = 1 TeV, X t = 0. Dashed blue lines provide rough 
estimates of the range of validity of the effective field theory. 



be seen in figure [2j We emphasize again that we actually do not need both coefficients c\ 
and C2 (corresponding to the U(l)y and SU(2)^ operators in (2.3)) to achieve enhancement, 
we could set e.g. C2 = 0, but the flexibility this additional parameter affords is useful in the 
next figure. The maximum allowed enhancement continues to decrease for larger values of 
tan /?, unless of course if we simultaneously lower M. 

From (3.17) we see that c^z^dec^ is maximized when c± and C2 have opposite signs. More- 



over, since is positive in (3.5), in order to achieve an enhancement of Rzj in (3.19), it is 
required that (02 — c\) > 0. This is seen in figure [2j where Rz^ is maximized for large positive 
values for C2 and large negative values for c±. Notice that the dependence on the sign of c\ 
and C2 for i? 77 and Rz-y is not specific to this scenario or SUSY. It simply follows from EW 
symmetry breaking, as can be seen in (3.6). 

In figures [3] and |4j we show a different representation of the same physics as in figure [2j 
where we fix the coefficients c\ and C2 in each curve, and instead vary the overall scale of 
new physics M. As expected, the effect of the higher-dimensional operators decreases with 
increasing M, but even for M approaching 10 TeV there can be some small effect. This 
perhaps somewhat counter-intuitive behavior is simply because the relevant SM couplings 
are small to begin with, as emphasized in the introduction. Since the "new physics" that 
generated these operators comes from a scale around or larger than 5 TeV, it will not be 
within easy reach of the LHC. 
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In figure [4] we have illustrated the behavior of R^z if we require for example i? 77 = 1.4 
(i.e. interpret the diphoton excess as signal) and vary — 5 < c\ < 5, so each curve represents 
a particular value of C2, and ends at some upper bound value of M where it is no longer 
possible to achieve the prescribed value of i? 77 = 1.4. We see that above tan (3 = 5 or so, one 
would have to rely on the scale M being not too far above 5 TeV. 




R 



M [TeV] 



tan (3 = 7 



c 2 = 1 



Co = -1 



M [TeV] 



Figure 4: The h — > rate Rz-y from (3.20) for various values of the coefficients 
ci and C2 of the operators in (2.2) and (2.3), if we require i? 77 = 1.4 and vary 
—5 < ci < 5 for fixed C2- The curve ends when c\ goes out of range for i? 77 = 1.4 
with the given parameters. We have set = 1 TeV, \i = 300 GeV, = 1 TeV, 
Xf = 0. Dashed blue lines provide rough estimates of the range of validity of the 
effective field theory expansion. 



With hindsight, it may appear that the analysis of the effects of O5 in (2.2) and Oq in 



(2.3) could have been performed mostly independently of each other. To be clear, we did 



not assume this: as a matter of principle, we always include the contribution to e.g. the 



mixing angle a from (2.2) when computing the effects of (2.3). But we emphasize that the 



i? 77 contributions arising from the dimension-six operators in (2.3) are maximized for small 



tan p, see (3.17). Therefore, since the usual MSSM tree level contribution to the Higgs mass is 



minimized for small tan /3, the contribution from the dimension-five operator to the tree level 
Higgs mass is crucial in order to accommodate a 126 GeV Higgs mass, as is seen in figure [T] 



5 Generating the effective operators from underlying physics 

The natural question is then what "new physics" could generate the effective operators dis- 
cussed. In this section we discuss a simple example of an underlying model from which both 
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05 an d O2 arise simultaneously in the low energy effective theory, upon integrating out some 
massive supersymmetric degrees of freedom. Consider a model that contains the following 
superpotential, involving a massive gauge singlet chiral superfielcj^] E: 



(5.1) 



and a gauge kinetic function r that depends on E, 



r(E) Tr (W a W c 



with r (E) D 



A 



(5.2) 



where A is a dimensionful suppression scale. If the SUSY mass /is is sufficiently large with 
respect to the energy scale under consideration, then E can be integrated out supersymmet- 
rically via its holomorphic equation of motion, which sets 



— H d ■ H u - -^-Tr {W a W a ) + 
MS MsA 



(5.3) 



where the dots stand for higher dimensional terms (further suppressed by fis, A). By inserting 
this solution back into the original Lagrangian, we obtain the following terms, 

A 2 



d z 6 fiH d • H u 



2/i S 



{H d ■ H u ) 



A 



H d -H u Ti(W a W a : 



M 


A 




Ms 



{H d -H u )\H d -H u ) 



(5.4) 



where we have included operators up to dimension-six. We see that operators O5, Oq of 



(2.2), (2.3) were simultaneously generated as a consequence of integrating out E. 



The dimension-six Kahler potential operator in the second line of (5.4) gives corrections 



to the quartic Higgs scalar potential, and hence to the tree level Higgs mass. However, in 



comparison to our dimension- five operator in the first line of (5.4), this operator is sup- 
pressed by a higher power of /x^. As long as fis is sufficiently large, the corrections from 
this dimension-six operator will be smaller in size with respect to the corrections from the 
dimension-five operator. 

It should be acknowledged that this example is not renormalizable since the gauge kinetic 
term (5.2) has dimension d = 5. To have a renormalizable microscopic model, one should also 
specify the degrees of freedom responsible for generating this d = 5 operator. Nevertheless 
this operator with a moduli-dependent gauge kinetic function is generically present in models 
derived from supergravity or string theory. 

Finally in order to connect with the discussion in the rest of the paper, we should assume 



that the dimension-five operators in (5.2) only involve the U(l)y and SU(2)^ gauge field 
strength, and not also the SU(3)c one that could in principle be also present. In string 
models, something along these lines could be achieved by for example considering a brane 



3 Gauge singlet fields with a supersymmetric mass term appear in general versions of the NMSSM [13II14IIT2 
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model in which the dimensionality of the branes that give rise to the U(l)y and SU(2) l gauge 
fields is different from those that give rise to the SU(3)c gauge field. In this way, since the 
U(l)y and SU(2)i branes and the SU(3)c branes, respectively, wrap different cycles of the 
internal geometry, they depend on different gauge singlet moduli fields associated with the 
different cycles. 

6 Conclusions 

Recent LHC data on the Higgs mass and its couplings to the photon, and the negative SUSY 
searches, present increasing difficulties for MSSM-like models to naturally accommodate a 
Higgs mass near 126 GeV, without undue fine-tuning, and a potential enhancement of the 
h — > 77 partial decay rate, without affecting other partial decay widths. Motivated by these 
observations, in this work we investigated whether supersymmetric effects beyond the MSSM 
(and even beyond reach for the LHC), could simultaneously accommodate these results. Using 
an effective approach, we identified two effective operators of dimensions d = 5 and d = 6 
that can address these problems and give the leading order corrections to the Higgs quartic 
coupling and the Higgs coupling to photons, respectively. 

We showed that the MSSM with small, supersymmetric corrections due to these effective 
operators can simultaneously naturally accommodate a Higgs boson with a mass near 126 
GeV, an enhanced Higgs coupling to photons (and also Z-y) relative to the SM expectation, 
and finally with SM-like Higgs couplings to the other SM particles. The scale of the super- 
symmetric effective operators is in the region of 5 TeV or even larger, and therefore possibly 
not within the LHC reach. The corrections from the dimension-six operators to the Higgs 
coupling to photons (and Zj) are maximized for small tan ft which is also the region where 
the dimension-five operator produces the most dramatic effect relative to the MSSM. This 
suggests that it is natural to consider these operators together and this is further supported 
by the fact that both of them can be generated simultaneously by an underlying model, as 
we showed. 

There remains the question of how the existence of these operators can be tested. Let 
us assume that the signal rate in the h — > 77 channel is confirmed to be higher than the 
SM expectation while the signal rates in all the other channels coincide with the SM values. 
If at the same time, one can rule out light stau sleptons in the mass range (of 150 GeV or 
so), needed in order to enhance the diphoton signal with the correct amount, this would 
cause a real problem for the MSSM and physics beyond the MSSM will be required. Should 
the excess go away when further data is analyzed, our results will remain useful to provide 
bounds on the scale M of supersymmetric "new physics" beyond the MSSM. Either way, this 
suggests that the diphoton rate is a useful, sensitive probe in this context. 
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A Appendix 

A.l Details concerning the Lagrangian 

Here we derive the Lagrangian of Section [2} Unlike in the text, we also include SUSY breaking 
effects associated with operators O5, Oq, by using the spurion field. The starting point is 

C =J d * e Yl (i-™^ 2 ) H U ViH i + {J ~d 2 0fi(l + B6 2 ) H d -H u + h.c.} +05 + 06 (A.l) 

i=u,d 

The superfield components are Vi = (Aj, Vi ifJ ,, Df/2), Hi = (hi,ipi, F.{). Also H d H u = e lJ H d Hi, 
with e« e k i = 8 ik ; e^e kl = 5 ik V l - 5 il V k , e 12 = 1, h d -h u = h° d h° u - h~h+. Further 

° 5= TiJ d2 ° (c ° + c '^ 2) {Hu ' Hd)2+ hx ' (A ' 2) 



= ^ [2 (h u -h d )(h u ■ F d + F u ■ h d -ip u ■ i>d)-(K ■ ip d +i[) u ■ h d ) 2 ] +^ (h u ■ h d ) 2 +h.c. 

and 

°e= ^2 E 1/2- / d2 ° ^ + c 's 92 ) MW a W a ) s (H u ■ H d ) + h.c. (A.3) 

s=1.2 9s K J 



A f p,vpa 

1 pa pa n 

tiv' 2 spv r spai 



s=l,2 

\/2 • + ^ • M(W + ^A^V) + (h u .F d + F u -h d -^ u - tp d ) X a s X a s } 



+ i^2(^-M(AX)+ h.c. 



Above we introduced T>^A a = \ a — g t ahc \ c for covariant derivatives of the gauginos. 
From C one finds the equations of motion for the auxiliary fields of Higgs superfields: 



p*i = _ e iP hP 



K- = -e«h' d [n + 2^ (h d -K)-^p^-^p\T\ (A.4) 
where q is a SU(2)^ doublet index. For the auxiliary fields of the vector superfields we find: 
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D a 2 



92 



{h\T a h d + hi T a h u ) (l - ( h u ■ hd + h.c.)) (A.5) 



C-2 



AM 2 



(h u -ijJd + tpu- h d ) A2 + h.c.) 



9i 



( h d(-\) h d + hl Q) (l 
V2 Cl 



ci 



4M 2 

with T a = cr a /2. The squares become 

C2 



2M 2 

(/i u • ip d + ipu ■ h d ) A" + h.c.) 



(K ■ h d + h.c/ 



D a 2 D a 2 



9l 



A 

V2 



1 -{w* hu ' hd + h - c - 

■g 2 \h d T a h d + hiT a h u 



\h d \ 2 - \h u \ 2 ) 2 + A\h\h^ 2 



_C2_ 

2M 2 



{h u -^d + ^u- h d ) \ a 2 + h.c. (A.7) 



(A.6) 



9l 
4 

V2 



gi 
2M 2 



h u -h d + h.c. 



-.91 



M - 

2M 2 



-l^l 2 ) 2 

{h u -ipd + ipu- h d ) Ai + h.c. 



(A.8) 



O5 and ©6 and eqs.(A.4) to (A.8) give the corrections to the MSSM Higgs Lagrangian. 



Using the corrected auxiliary fields in the usual MSSM Higgs Lagrangian, additional terms 
suppressed by 1/M and 1/M 2 are generated. The full on-shell Lagrangian is then: 



£ = £d + Cf + Ci + £-2 + £3 + -Cssb • 



Eliminating the D-dependent terms in £ one finds, see eqs.(A.5) to (A.8): 



£ D = Y,-l D s D s[ l + ^p( h -- h <i + h - C -)] 



8=1,2 



(A.c 



(A.10) 



and use (A.8). Eliminating the -F-dependent terms in £ gives £f~. 



£ F = \F d \ 2 + \F u \ 2 = \fi + 2^h d -h u \ 2 {\h d \ 2 + \h u \ 2 ) 



+ 



^{-■m^-m^ (N 2 + i/g 2 )+h.c. 



(A.ll) 



Apart from auxiliary field contributions, there are also terms in the Lagrangian with space- 
time derivatives, that contribute to the kinetic terms for Weyl fermions ip u di Af 2 when the 
neutral singlet h u d components of h Ut d acquire a vev: 



£s 



C2 

AM 2 



{h u ■ h d ) [i {\ a 2 a»V u \ a 2 - P M A^A§)] + h.c. +(2^1) 



(A.12) 
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When the Higgs neutral singlets acquire a vev, these terms produce wavefunction renormal- 
ization of Weyl kinetic terms and a threshold correction to gauge couplings g\ and g 2 . 

There are also terms contributing to fermion masses when the Higgs fields acquire vev's 



= ^ 2 (h u -h d ){X a 2 X a 2 ) + ^{h u -h d ){X 1 X l ) 

+ jjji[2 (h u ■ h d )(-if) u ■ ip d ) - (h u -ipd + ^u- h d ) 2 ] + h.c. 



(A.13) 



Further, there are some interaction terms generated 



_C2_ 

AM 2 



\ (K ■ h d ) (F^F^ V + \ ^F 2 ;^ pff; 



V2 {h u -iPd + ^u- h d ) <J» v X a 2 F^ v -^u-^d AgAg] + (2 -> 1) + h.c.} (A.14) 



Finally, the Lagrangian contains {F and D independent) corrections from supersymmetry 
breaking due to spurion dependence in the dimension 5 operator as well as the usual soft 
terms of the MSSM. All these together give Cssb'- 



£ssb = -Vssb 



M 



{h u • h d f + LiB{h d - h u )+ h.c. - m 2 d \h d \ 2 -m 2 u \h u 



This concludes the presentation of the Lagrangian to 1 /M 2 order. From C we find the scalar 



potential Vh for the Higgs sector shown in the text, eq. (2.6) in which as usual h u>d denote 
SU(2) doublets. From this one obtains: 



m i 



loop + Aml 



(A.15) 



with: 



nr 



h, loop 



■»' 



1 f O O i — "1 

- yn A + m z + 5 m z sin /3 — y/w J 



[{m\-m z ) cos 2/3 + 5 mf sin^ /3f + sin^ 2/3 {m\ + m z y (A. 16) 



and 



Am 2 = f l {2^) + f 2 {-2 C ^) + h{2^) 2 + h{-2 C ^) 2 + h{2^){-2 C ^) 



(A.17) 



where 
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h 
h 
h 

h 
h 
h 



v 2 sin 



2/3 { 



1 + 



(m 2 + m|) 



} 



cos 2/3 



w 



' sin 2 2/3 + 4 



4 M 



(m\ — m? z ) cos 2/3 + m? z 5 sin 2 /3 } 

2/3} 



| - 1 + _( m 2 +m | J sjl p 



+ 



W 



3/2 



(m^ + m|) 2 v A sin 2 2/3 



(777/4 + m|) 2 sin^4/3 



16 u> 3 / 2 



32 



4w 3 / 2 
sin 2/3 



(m^ + m|) [ S m 2 z + (2 777/4 - (2 + 5) m|) cos 2/3] sin 4/3 



1 + 



1 



[8m 2 4 -(4+3(5)m|+6 ( 5m|cos2/3+3(4m 2 1 - <5m|) cos4/3]l (A.18) 



and finally 

2 2t5/j 
m ^ = s4n^ 



4 2 

7J COS 



2^ (sfo £25|A +0 /J_ ' 
32 sin2/3 \ M 2 M 2 / V M3 , 



A. 2 1-loop form factors 



The form factors in (3.4) are given by, 



Af = \A 1/2 {r t ) , Af = \A l/2 {n), 
A™ = At{rw) , Af = N c Q 2 A l/2 (r t ) , 

Af^^A^Aw) , A% = N c 9^t ~ 4 -9^M Al/2 (r t ,Xt) (A.19) 

COS U yj 

where n = 4m 2 /m 2 h , X { = 4m 2 /m 2 z , N c = 3, Q t = 2/3, T 3 W = 1/2 and 

A 1/2 (t) = 2T 2 [r- 1 + (T- 1 -l)f(T- 1 )] , 
Ai(t) = -r 2 [2r- 2 + 3T- 1 + 3(2r- 1 -l)/(T- 1 )] , 

A 1/2 (T,X) = h(T,\)-I 2 (T,X) , 

A^t, A) = 4(3 - tan 2 9 w )I 2 (t, A) + [(1 + 2t~ 1 ) tan 2 6 W - (5 + 2t~ 1 )] I x {t, \\k.2Q) 

where 

/i(r,A) = ^ + ^r)-rt^r I )-j(Al 1 

^2(r,A) = --^^[/(r-i)-/^ 1 )] , (A.21) 
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and 



/(*) 



g( x ) 



arcsm Wx 



log 



1 + VT 



27T 



yV 1 - 1 



arcsin \/x 



log 



l + VT 



„-i 



X < 1 
X > 1 , 



(A.22) 



1-vT 



„-i 



ITT 



X < 1 
X > 1 



(A.23) 
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